Noise and Counting Statistics of Insulating Phases in One-Dimensional Optical 

Lattices 
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We discuss the correlation properties of current carrying states of one-dimensional insulators, 
which could be realized by applying an impulse to atoms loaded onto an optical lattice. While the 
equilibrium noise has a gapped spectrum, the quantum uncertainty encoded in the amplitudes for the 
Zener process gives a zero frequency contribution out of equilibrium. We derive a general expression 
for the generating function of the full counting statistics and find that the particle transport obeys 
binomial statistics with doubled charge, resulting in super-Poissonian noise that originates from the 
bosonic nature of particle-hole pairs. 



The characterization of particle fluxes using noise and 
full counting statistics (FCS) has played a very prominent 
role in quantum optics and, more recently, mesoscopic 
physics. Lately, these same ideas have come to the fore in 
the study of cold atomic gases, where noise measurements 
have been shown to be sensitive probes of inter-particle 
correlations 0, Q as well as quantum statistics 0, OJ ■ 

Although the idea of measuring the full distribution 
of particle transport has generated great theoretical in- 
terest recently for the case of massive particles (usually 
electrons) 0.[g[. the experimental determination in Ref.y 
of the FCS for an atom laser represents the first such 
measurement. The observation of Poissonian statistics 
without bunching confirms that the atom laser shares 
many characteristics with its optical counterpart. This 
means that the noise spectrum 



S(u)= / dte iwt (J(0)J(t) + J(t)J(0)>, 



(1) 



is proportional to the current: S(u>) = 2(J). The zero 
frequency contribution is usually called shot noise, and 
is characterized by a Fano factor F = 5(0)/2(J), here 
equal to 1. The FCS is characterized by the generat- 
ing function x(A, T) for the moments of the number of 
particles transported in time T. Its logarithm generates 
the cumulants of the distribution. Thus for the Poisson 
distribution we have 

lnx(A,T) = (J)T(e lA -l). 

Most discussions in mesoscopic physics concern the gen- 
eration of shot noise from partitioning of the particle 
flux by localized scatterers. In this Letter we consider 
another natural source of shot noise: the Zener pro- 
cess by which current-carrying states are generated on 
a lattice [jj. Zener tunneling has been observed in one- 
dimensional optical lattices loaded with cold atoms ||, 
and this is a natural setting for the phenomena that 
we will discuss. Our aim is to extend Zener's result 
— ($/27i7i)m[l — Plz] for the transition rate per unit 
volume, where Plz is the Landau-Zener probability and 
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FIG. 1: a) Schematic view of particle-hole pair production 
in a one- dimensional insulator subject to a force, b) The 
discretization used to evaluate the counting statistics. 



$ the applied force, to the entire distribution of the re- 
sulting charge transport. 

When fermions are partitioned by a localized scatterer 
with transmission probability T, it was found that the 
Fano factor is 1 — T :9.j, i.e. the noise is suppressed rel- 
ative to the Poisson case. In Ref. this suppression 
was revealed to be a consequence of the binomial statis- 
tics of charge transport. In the present case, we might 
expect by analogy that when Plz is small, Poissonian 
statistics is obtained with a Fano factor of one, falling 
to zero as the probability rises to unity. Such an expec- 
tation corresponds to 'classical' picture of particle-hole 
creation, illustrated in Fig. QJi). Remarkably, this ex- 
pectation is only partially borne out. We find that the 
FCS resulting from an impulse applied to an insulating 
state on a one-dimensional lattice is binomial, but with 
doubled 'charge'. This means that the Fano factor is 
F = 2 (1 — Plz)- The origin of this unusual 'bunching' 
behaviour, which can occur for bosons or fermions, is the 
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coherent creation of particles and holes, or alternatively 
the bosonic behavior of particle-hole pairs. 

There are two natural situations to consider. Fcrmions 
loaded onto a lattice form a band insulator at integer fill- 
ing. In the absence of interactions current-carrying states 
can be long-lived as there is no way for quasi-momentum 
to relax [l(|. Restricting ourselves to two bands, one 
empty and one full - more bands present no difficulty - 
we fix for concreteness the massive Dirac band structure 
near the Brillouin zone boundary 



Cp 03 + U Oi . 



(2) 



Here p is the deviation of the (quasi-)momentum from 
the Brillouin zone boundary, c is a velocity, set to unity 
from now on, and 2u is the bandgap, proportional to the 
intensity of the laser light field that mixes left and right 
moving particles. In the basis of left and right movers, 
the two bands have the form 



l/)-(p) 



1 



^2E p {E p -p) \ E p -p 

1 (p-E, 
y/2E p (E p -p) 



(3) 



with energies E±(p) = ±E p = zty^p 2 + u 2 . Introducing 
the corresponding creation and annihilation operators, 
the normal-ordered current operator at the origin can be 
expressed as 



J(o) = J2 v +(p>p'^Up)Mp') 

Ptp' 

-v B (p,p')^_(p')^Up) 

+ \w B {p,p')^ + {p)4,_{p') 



(4) 



(we work in units with H = 1) Here 

u 2 ~{E p -p) (E p , -p') 



v B (p,p') = 
v B (p,p') = 
w B (p,p) = 



2^E p (E p ~p)^E v AE p 



■p>) 



-v?(p,p') 



u(E p -p + E p , - p') 



2^E p {E p -p)^E p ,(E p , -p>) 



(5) 



When p = p' the band diagonal parts of the cur- 
rent Eq. involve the group velocities v B (p,p) = 
dE ± (p)/dp. 

We will also consider a one-dimensional bosonic Mott 
insulator, modeled by hardcore bosons b 2 — on a lattice 
at half filling, with a staggered field 



H 



M 



h.c 



(G) 



This model is equivalent under a Jordan- Wigner trans- 
formation to a band insulator of fermions with two bands. 



The expression for the current, found by standard meth- 
ods has the same form as Eq. (0J, with 



(E p , -u)£ p e* +(E p -u)^ p ,e- i P 



2y/E p (E p -u)^E p ,(E p , -u) 
v M (p,p') = -v?(p,p') 



w M (p,p') = 



^e'P ~{E p ~u){E p , -u), 



2y/E p (E p -u)y/E p .(E p . -u) 



; (7) 



with E p = y/£ p + u 2 and £ p = sinp, with p the deviation 
from 7r/2, the zone boundary. For p< 1 Eq. JJJ coincides 
with Eq. J^}. The specific forms of v±(p,p') and w(p,p') 
will not be needed. We now consider what happens 
when an impulse is applied to the lattice, that is, a time- 
dependent but spatially constant force <&(<). The prob- 
lem remains translationally invariant with each momen- 
tum state Pi being mapped to one at pf — Pi + J ^(t')dt' . 
In general, there is some amplitude to make a transition 
to the upper band, meaning that the wavefunction of a 
given momentum state, in the basis of upper and lower 
bands, evolves as 



*/(P/) 



B, 



A 



We now use these amplitudes to compute the resulting 
current fluctuations. Let us begin by finding the aver- 
age current. Taking the expectation of the Heisenberg 
representation form of J(0, t) 



{J(0,t))=J2[v+(P,P)-V-(j>,p)]\B p 



+2w{p,p)ReB;A p e l( - E +^~ E - 



(p))t 



(8) 



The last term gives rise to a transient, and in the t — > oo 
limit we have the natural result 



(J(0,t 



dE+jp) dE_(p) 
dp dp 



\B„ 



(9) 



Moving on to the current fluctuations, the expectation 
for the connected correlation function 

(j((M)j(o,t'))c = (J(o,t)j(o,*')> - (J(o,t))(J(o,t')), 

is straightforwardly obtained but rather lengthy. It is 
useful to first consider the equilibrium case when A p — 1 , 
Bp — 0. Then we have 



-i(E+(p)-E-(p'))t 



(j(o,*)j(o,o)) c = j/Mp',p)< 

p,p' 

(10) 

Note that this arises from off-diagonal part of the cur- 
rent Eq. Q. For the Dirac case the integrals are readily 
performed to give 

(J(0,<)J(0,0)) c = ^ [K^-iutf 

+ (K (-iut) + TrH (ut)) 2 ] ,(11) 
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where K n (x) is a Bessel function and H n {x) a Struve 
function. At t — > we have the behavior 



(J(0,t)J(0,0)) c 



1 



2TT 2 (t-ioy 



In this limit the correlator matches that of the massless 
system, which is fixed by the Schwinger commutators 



Jl/r{%),Jl/r{%') 
Jl/r(x),Jr/l(x') 



0. 



(12) 



As may be readily seen from Eq. QlOfl. the noise spec- 
trum contains only positive frequencies, indicating that 
the system can only absorb energy at zero temperature. 
Furthermore the spectrum is gapped at 2u, meaning that 
there is no noise at zero frequency. 

It is thus reasonable to drop the gapped contribu- 
tions to the noise (we will justify this more rigorously 
later). Additionally, we now focus on the particle- hole 
symmetric case where E±(p) = iE p , and v+(p,p') — 
—V—(p,p?) = v(p,p'), as in our two examples. Then the 
current correlator out of equilibrium is 



(J(0,t)J(0,t')) c = 2j2 »(p,p>(p^)[e-^*-^)(*-0^ (13) 
p,p' 



Notice that the noise correlator is non-stationary, since 
the non-equilibrium state is not an eigenstate. It is still 
possible to define a Fano factor as {Q%) c = F(Qt), where 
Qt is the current passed in time T. Then Eq. JJJJ and 
Eq. Ijl3|l give (we use the property v(p, —p) = 0) 

2E P 



In the above we have introduced the doublet 



F = 



\E'{p)\\A p \ 2 \B^ 



(14) 



E p E>(p)\B p \i 

This result closely resembles the formula for the shot 
noise of a point contact. The key difference lies in 
the factor of two, so that in the limit B p — ► we get 
F = 2 rather than 1. The origin of this surprising super- 
Poissonian behavior is the second term of Eq. (|13J) . 

We now discuss the FCS for this problem. Two toy 
models of 'ideal' detectors have been considered in the 
literature [l2l - both give the generating function of 
the FCS as 



X(A,T) 




dt'sgn c (t')J(0,t')) ).(15) 



Here T denotes the measurement time, and c(T) is a 
(Keldysh) time contour that passes from to T and back 
again. T c is the operation of time ordering on this con- 
tour, and sgn c is +1 on the forward branch, and — 1 on 
the backward branch. The expectation value is taken 
over the current-carrying state of the system. It should 
be noted that the measurement performed in Ref.0 would 
seem to have little in common with the idealized detec- 
tion arrangements for which the above form of x(A, F) 
was found. We hope to provide a detailed analysis of 
this issue in a forthcoming publication: for the moment 
we will work with Eq. I|15l) for want of something better. 

We use Wick's theorem to rewrite Eq. I|15|) in normal- 
ordered form 



X (A,T) =xvac(A,T)( :exp(- 
XVac(A,T) = exp [Trln^r 1 ] , 



-zTr 



G = (t 



> 



iAV)~\l6) 



4>-(p) 



while 2J PiP ,(t,f) = T c * p (<) ® *J,(f)- : * p (t) ® ¥ pl (t') 
is the Feynman propagator 



vpotw) °p,p> y o -e e (t',t) J 



9 c (t,t') is 1 if t comes after t' on the contour, and zero 
otherwise. The matrix A p y(t) is introduced so that 



A 



sgn c (t)j(0,t) =J2H A P.p'(t)* P >- 



That is, it contains the v±(p,p') and w(p,p') of Eq. J3J. 
The 'Tr' in Eq. (|16|l is a trace over the band, time, and 
p indices. 

In finding Eq. (|16fl we have factored off the part of the 
generating function that exists in the 'vacuum'. Recall 
that we earlier showed that the quantum noise Eq. 1|10|) 
in the ground state has a gapped spectrum. We thus 
do not expect it to contribute to the FCS at long times 
T. Since Eq. I|1U|) originated from the band off-diagonal 
contribution to the current, we can drop these effects at 
the outset by setting w(p,p') = in A PiP /. 

Alternatively, we can perform a clean calculation in 
the following way ^3] • We pass to a set of equally spaced 
states in energy space with energies ±E n = ±it ± 2im/F 
(assuming also that 2u is an integer multiple of 2w/T). 
In this basis the current operator takes the form 



4 



j(o,t)= 



v bb >(n,n r ) 



,b,b> 



\vbb(n, n)v b i b >(ri, n')| 1/2 



V+bn'y+b'n' t ^ V-b'n'V-bn^ 



(17) 



The 'branch' index & = — 1, 1 is needed to account for p < and p > states with momenta bp n such that E n — E(p n ) 
(see Fig.^.). In Eq. (|17|l v bb i(n,n') = v(bp n) b' p n >). The advantage of this approach is that not only do the effects 
of w(p,p') disappear but x(A,T) takes a diagonal form with Xvac (A,T) = 1 and 



x(A,T) = (:exp (t£*L 



7 J , 6 











bn 



(18) 



Eq. I|18l) represents a general form for x(A, T), ready for averaging over the state of the system. The state following 
an applied impulse is 



H[Abn + B b Jl b J_ bn ] |0>. 



(note that the assumption of particle-hole symmetry enters here in an essential way - pairing of particle and hole 
states at fixed momentum translates to pairing at given n) This state has 



— bn/ 
bn / 



l 



T \ \ n ^bn 

The average of the normal-ordered exponent in Eq. I|18|) can then be found 



\B bn \ 2 



(19) 



X (A,T) = JJdet 

1 1 , b 

ln X (A,T) = Y,hx[\A 



2 + \B bn \ 2 e* xb B^A 



iXb 



bn bn 



I Am 

AbAil-e^) \A b7 f + \B bn 



2giAfc 



We used the normalization condition |^4 bra | 2 + \B bn \ 2 = 1. 
To make things concrete, let's consider the case of a con- 
stant force $ applied for time r, with the resulting im- 
pulse $r not creating particle-hole pairs too far away 
from the zone boundary, so that the Dirac structure 
Eq. J2J is a good approximation. As long as $r s> u/c, 
however, most of the resulting particle-hole excitations 
are located on the linear part of the spectrum at p > 0, 
so that we can neglect the crossing point and take the 
amplitudes to be given by the Landau-Zener formula 

2 JO b = -l 

\ B bn\ - < 

yp LZ = e <=* = 1, 

giving the FCS 

X (A,T)= [(l-P LZ )+P LZ e^] c * rT/2 \ (20) 

and the Fano factor F = 2 (1 — Plz), in agreement with 
Eq. (fTljl . The result Eq. f2T)| corresponds to binomial 
statistics of particle transport, with a probability Plz of 
transporting two charges (a particle passing in one direc- 
tion and a hole in the other) per 'attempt', and a number 



of attempts cQtT/2-k. Note that the unusual doubling is 
a consequence of the coherent creation of particles and 
holes: if the off-diagonal parts of Eq. (|19f) are set to zero, 
the FCS is a product of binomial statistics for particles 
and holes separately, with Fano factor < 1 . The classical 
picture of Fig QJi) upon which our intuition was based 
corresponds to this case. 

A brief comment is in order concerning the case with- 
out particle-hole symmetry. In general the second term 
in Eq. i|13|l will give rise to a zero frequency contri- 
bution to the noise from a pair of momenta p,p' with 
E + (p) + E-(p') — 0. This contribution involves the 
combination of amplitudes B*A p B p ,A*,, and is thus in 
general sensitive to the phase (the 'Stokes phase'). One 
would expect a general decrease of the Fano factor for this 
reason, though super-Poissonian behavior is still possible. 

To conclude, we have investigated the counting statis- 
tics of the Zener process, and discovered an unusual 
super-Poissonian behavior. It would be interesting to 
perform such calculations for more realistic models of in- 
sulators, where particles and holes can interact, to see 
how the results are modified. 
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